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Abstract 
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Ground-state properties of ^ = 20 nuclei (^"N, ^'^O, ^"F, ^"Ne, 



j23 . ^°Na, ^^Mg) are described in the framework of Relativistic Hartree- 

Bogoliubov (RHB) theory. The model uses the NL3 effective interac- 
tion in the mean-field Lagrangian, and describes pairing correlations 
by the pairing part of the finite range Gogny interaction DIS. Binding 

p^ \ energies, quadrupole deformations, nuclear matter radii, and differ- 

ences in radii of proton and neutron distributions are compared with 
recent experimental data. 



PACS: 21.60.Jz, 21.10.Dr, 21.10.Gv, 27.30.+t 



1 Introduction 

Sizes, shapes and binding energies are fundamental characteristics of nuclei 
and reflect the basic properties of effective nuclear forces. The description 
of the generalized moments of the nuclear density distributions provides an 
important test for nuclear structure models. In particular, a description of 
ground-state properties of an isobaric sequence of nuclei tests the isovec- 
tor channel of the effective nuclear force. The correct parameterization of 
this channel is essential for the description of structure phenomena in exotic 
nuclei far from /^-stability. In the present study we analyze the sequence 
of A = 20 nuclei: ^on, 209, 20p^ 20^^^ 20Na, ^OMg. Experimental data on 
nuclear matter radii derived from the measured interaction cross section of 
these nuclei have been reported in Ref. |l|]. Evidence has been found for the 
existence of a proton skin for ^°Mg and of a neutron skin for ^°N. The largest 
difference in radii (~ 0.2 fm) has been reported for the mirror nuclei ^°0 and 
2°Mg. This last result has prompted several theoretical investigations, but 
the extremely small value for the matter radius of ^°0 has not been repro- 
duced in calculations which include the microscopic three-cluster model [Q], 
the Hartree-Fock model with Skyrme interactions 0, ^, ^, |^ , a single parti- 
cle potential model 0, and the shell- model 0. The experimental data on 
interaction cross sections of A = 20 isobars were also reviewed in Ref. [0] by 
using a variety of phenomenological and semi-microscopic models. 

In the present analysis of the A = 20 isobaric sequence we employ the Rel- 
ativistic Hartree-Bogoliubov (RHB) theory. Based on the relativistic mean- 
field theory and on the Hartree-Fock-Bogoliubov framework, the RHB theory 
provides a unified description of mean-field and pairing correlations. It has 
been successfully applied in the description of nuclear structure phenomena. 



not only in nuclei along the valley of /^-stability, but also of exotic nuclei 
close to the particle drip lines. In particular, applications relevant to the 
present study include: the halo phenomenon in light nuclei ||], properties of 
light nuclei near the neutron-drip 0, reduction of the spin-orbit potential 
in nuclei with extreme isospin values |T^, the deformation and shape coex- 



istence phenomena that result from the suppression of the spherical N=28 
shell gap in neutron-rich nuclei [|ll|], properties of proton-rich nuclei and the 



phenomenon of ground-state proton radioactivity |12, 13 



nu- 



The relativistic mean field theory is based on simple concepts [0 
cleons are described as point particles, the theory is fully Lorentz invariant, 
the nucleons move independently in mean fields which originate from the 
nucleon-nucleon interaction. Conditions of causality and Lorentz invariance 
impose that the interaction is mediated by the exchange of point-like effective 
mesons, which couple to the nucleons at local vertices. The single-nucleon 
dynamics is described by the Dirac equation 



,0 , „ ^ J) , ^(1 ~ ^3) .0 



<^ -ioL ■ V + (3{m + g^a) + g^uj^ + gpT-^pl + e A^ \ ipi = e^ipi. (1) 

0", uj, and p are the meson fields, and A denotes the electromagnetic potential. 
(?o- gui, and gp are the corresponding coupling constants for the mesons to the 
nucleon. The lowest order of the quantum field theory is the mean-field 
approximation: the meson field operators are replaced by their expectation 
values. The sources of the meson fields are defined by the nucleon densities 
and currents. The ground state of a nucleus is described by the stationary 
self-consistent solution of the coupled system of Dirac and Klein-Gordon 
equations. 

In addition to the self-consistent mean-field potential, pairing correlations 
have to be included in order to describe ground-state properties of open-shell 



nuclei. In the framework of the relativistic Hartree-Bogohubov model, the 
ground state of a nucleus |$ > is represented by the product of independent 
single-quasiparticle states. These states are eigenvectors of the generahzed 
single-nucleon Hamiltonian which contains two average potentials: the self- 
consistent mean-field T which encloses all the long range particle-hole (ph) 
correlations, and a pairing field A which sums up the particle-particle {pp) 
correlations. In the Hartree approximation for the self-consistent mean field, 
the relativistic Hartree-Bogoliubov equations read 

fhn-m-\ A \fU,{v)\_fU,ir)\ 

[ -A* -h^ + m + XjWkir)) ""'Wkir))- ^^^ 

where ho is the single-nucleon Dirac Hamiltonian (0), and m is the nucleon 
mass. The chemical potential A has to be determined by the particle num- 
ber subsidiary condition in order that the expectation value of the particle 
number operator in the ground state equals the number of nucleons. A is the 
pairing field. The column vectors denote the quasi-particle spinors and Ek are 
the quasi-particle energies. The RHB equations are solved self-consistently, 
with potentials determined in the mean-field approximation from solutions 
of Klein-Gordon equations for the sigma meson, omega meson, rho meson 
and photon field, respectively. 

2 Ground-state properties of A=20 isobars 

The details of the ground-state properties of the A = 20 isobaric sequence will 
depend on the choice of the effective nuclear interaction. In the last few years 
several new and more accurate parameter sets of meson masses and meson- 
nucleon coupling constants for the effective mean-field Lagrangian have been 
derived. Due to large uncertainties in the experimental data, older effective 
interactions were not specifically designed to describe ground-state isovector 



properties. Only more recently have isovector properties been included in the 
set of data on which the effective interactions are adjusted. In the relativistic 
mean-field model, perhaps the most accurate set of meson masses and meson 
nucleon coupling constants is NL3 |TB]. In most recent applications of the 



RHB theory we have used the NL3 effective interaction for the mean-field 
Lagrangian, and the pairing field has been described by the pairing part of 
the Gogny interaction with the DIS parameter set |16|. This force has been 



very carefully adjusted to the pairing properties of finite nuclei all over the 
periodic table. In particular, the basic advantage of the Gogny force is the 
finite range, which automatically guarantees a proper cut-off in momentum 
space. 

The RHB equations are solved self-consistently, with potentials deter- 
mined in the mean-field approximation from solutions of Klein-Gordon equa- 
tions for the meson fields. The Dirac-Hartree-Bogoliubov equations and the 
equations for the meson fields are solved by expanding the nucleon spinors 
Uk{r) and Vfc(r), and the meson fields in terms of the eigenf unctions of a 
deformed axially symmetric oscillator potential. The calculations for the 
present analysis have been performed by an expansion in 12 oscillator shells 
for the fermion fields, and 20 shells for the boson fields. A simple block- 
ing procedure is used in the calculation of odd-proton and/or odd-neutron 
systems. The blocking calculations are performed without breaking the 
time-reversal symmetry. A detailed description of the Relativistic Hartree- 



Bogoliubov model for deformed nuclei can be found in Ref. [|13 |. 

In Fig. |l] we compare the RHB NL3+D1S binding energies of the A = 20 
nuclei (^^N, ^OQ, ^op, ^ONe, ^^Na, ^"Mg) with experimental data [0- The 



agreement between theory and experiment is very good. The self-consistent 



RHB neutron and proton ground-state density distributions for the A = 20 
isobaric sequence are shown in Fig. |^. We will eventually compare the cal- 
culated radii with experimental values, but here we notice the trend: a pro- 
nounced proton skin is observed in ^°Mg, it slowly disappears and the proton 
and neutron density distributions are almost identical in ^''Ne, the neutron 
skin develops and is very pronounced in ^°N. This trend is in complete agree- 
ment with experimental evidence for the existence of a proton skin for ^°Mg 
and of a neutron skin for ^°N [|l|. 

The calculated ground-state quadrupole deformations of A = 20 isobars 
are shown in Fig. ^ as function of the isospin projection T^. Except for 
^°Ne and ^''Na, the A = 20 nuclei are essentially spherical. ^°Ne is strongly 
prolate deformed {P2 = 0.365), and an oblate ground-state is calculated for 
^°Na {P2 = —0.153). The surface thickness and diffuseness parameters of the 
proton and neutron density distributions (Fig. |l|) are displayed in Fig. ^ The 
surface thickness s is defined to be the change in radius required to reduce 
p(r)/po from 0.9 to 0.1 (po is the density in the center of the nucleus). The 
diffuseness parameter a is determined by fitting the neutron density profiles 
to the Fermi distribution 

Pir) = Po f 1 + exp{ — ^) j ) (3) 

where Ro is the half-density radius. The surface thickness and diffuseness 
parameters quantify the proton and neutron skins: for ^°Mg s{p) = 1.81 fm, 
s{n) = 1.50 fm, and the corresponding diffuseness parameters a{p) = 0.42 
fm, a{n) = 0.35 fm; for ^ON s{p) = 1.60 fm, s{n) = 1.94 fm, a{p) = 0.38 fm, 
a{n) = 0.46 fm. For ^°Ne the surface thickness and diffuseness parameters 
of the proton and neutron distributions are identical. 

The calculated radii are compared with experimental data |]l], |l^ in Figs. ^ 
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and ^. In Fig. ^ we display the nuclear matter radii (left panel), and the 
differences between proton and neutron radii (right panel) of the A = 20 
nuclei as functions of the isospin projection T^. The theoretical values are in 
very good agreement with experimental data. The calculated matter radii 
reproduce the observed staggering between even and odd values of T^. The 
only serious discrepancy is ^°0: the calculated matter radius is 2.79 fm and 
the experimental value is 2.69(3) [|l|. However, as we have already mentioned 
in the introduction, several theoretical models have failed to reproduce this 
extremely small value for the matter radius of ^''O p, ^, 0, 0, and it was also 
suggested that the reported experimental value should be reconsidered [§]. 
All theoretical calculations predict the matter radius in the interval between 
2.77 fm and 2.83 fm, in agreement with the result of the present analysis. 
In fact, since most models predict over-binding for the O isotopes (^^O, and 
even ^^O, are bound nuclei in most model calculations), one would expect 
the theoretical value for the matter radius in ^''O to be smaller than the 
experimental one. With the exception of ^''O, also the calculated differences 
between proton and neutron radii (right panel) are in very good agreement 
with experimental data. In particular, the theoretical values agree with the 
large differences found in nuclei with a large excess of neutrons (^°N, 0.33 ± 
0.15 fm) or protons (^°Mg, 0.50 ± 0.28 fm). 

Finally, in Fig. ^ theoretical and experimental values are compared for 
the proton radii as function of the isospin projection T^ (left panel), and for 
the matter radii as function of the binding energy. We notice a very good 
agreement for all the proton radii. The RHB model also reproduces, with the 
above discussed exception of ^°0, the functional dependence of the matter 
radii on the binding energy of the A = 20 nuclei. 



3 Summary 

In this study we have apphed the Relativistic Hartree-Bogohubov theory in 
an analysis of ground-state properties of nuclei that belong to the A = 20 
isobaric sequence. The NL3 effective interaction has been used for the mean- 
field Lagrangian, and pairing correlations have been described by the pairing 
part of the finite range Gogny interaction DIS. This particular combination 
of effective forces in the ph and pp channels has been used in most of our 
recent applications of the RHB theory. Theoretical predictions for binding 
energies, neutron and proton ground-state density distributions, quadrupole 
deformations, nuclear matter radii, proton radii and differences between pro- 
ton and neutron radii for ^°N, ^°0, ^°F, ^°Ne, ^°Na, and ^°Mg have been 
analyzed and compared with available experimental data. It has been shown 
that the NL3 effective force provides a very good description of the observed 
ground-state properties as function of the isospin projection T^. The results 
of the present analysis confirm that the isovector channel of the NL3 inter- 
action is correctly parameterized and that this effective force can be used to 
describe properties of nuclei far from /^-stability. 
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Figure 1: Binding energies of A = 20 nuclei calculated with the NL3 + 
Gogny DIS effective interaction. The theoretical values are compared with 



the experimental binding energies |T7 |. 



Figure 2: Self-consistent RHB neutron and proton ground-state density 
distributions of A = 20 nuclei. 



Figure 3: Calculated ground-state quadrupole deformations of A = 20 
isobars as function of the isospin projection T^. 



Figure 4: Surface thickness of the neutron and proton density distributions 
of the y4 = 20 isobars, calculated with the NL3 + Gogny DIS effective 
interaction. In the insert the corresponding values for diffuseness parameter 
are included. 



Figure 5: The nuclear matter radii (left), and the differences between proton 
and neutron radii (right) of the A = 20 nuclei as functions of the isospin 
projection T^. Results of fully self-consistent RHB calculations are compared 
with experimental data [|I|. 



Figure 6: Theoretical values for the proton radii as function of the isospin 
projection T^ (left), and for the matter radii as function of the binding energy 
(right) of the A = 20 isobaric sequence, in comparison with experimental 
data M. 
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